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We study theoretically the properties of SIFS type Josephson junctions composed of two superconducting (S) 
electrodes separated by an insulating layer (I) and a ferromagnetic (F) film consisting of periodic magnetic 
domains structure with antiparallel magnetization directions in neighboring domains. The two-dimensional 
problem in the weak link area is solved analytically in the framework of the linearized quasiclassical Usadel 
equations. Based on this solution, the spatial distributions of the critical current density, Jc, in the domains 
and critical current, 7c, of SIFS structures are calculated as a function of domain wall parameters, as well as 
the thickness, dy, and the width, W, of the domains. We demonstrate that Icidy, W) dependencies exhibit 
damped oscillations with the ratio of the decay length, ^i, and oscillation period, ^2, being a function of the 
parameters of the domains, and this ratio may take any value from zero to unity. Thus, we propose a new 
physical mechanism that may explain the essential difference between and ^2 observed experimentally in 
various types of SFS Josephson junctions. 


It is well known that properties of Josephson struc¬ 
tures with ferromagnetic (F) material in a weak link 
region depends on relation between the complex de¬ 
cay length, -I- i^2^) geometrical 

parameters of these junctions n-0- If F metal is in 
the dirty limit and exchange energy, 77, sufficiently ex¬ 
ceeds the critical temperature of superconducting (S) 
electrodes, ttTc, then from Usadel equations it follows 
that ~ ^2- However, it was demonstrated experimen¬ 
tally m-m that there could be a noticeable difference 
between and ^2- Previously the difference has been 
attributed either to the presence of strong paramag¬ 
netic scattering in the F layer [7], or to violation of 
the dirty limit conditions in ferromagnetic material 
[T^ . However, application of the first of the mechanisms 
for the experimental data interpretation requires the ex¬ 
istence of unreasonably strong paramagnetic scattering 
in the weak link material [7]. The relation between an 
electron mean free, £, and ,^1, ^2 in typical experimen¬ 
tal situation is also closer to the dirty limit conditions, 
7 < ^1, ^2 rather than to the clean one. 

In this article we prove that the existence of a fer¬ 
romagnetic domain walls in F layer can also lead to 
appearance of substantial differences between ^1 and 
^2 even in the absence of strong scattering by paramag¬ 
netic impurities, and under the fulfilment of the dirty 
limit conditions in the F material. 
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FiglU Geometry of the considered SIFS Josephson 
junction and its enlarged part, which includes two 
halves of domains and domain wall separating them. 
The insulating barrier I has a small transparency 
(shown by a blue line). 

Model. Consider multilayered SIFS structure pre¬ 
sented in Fig.I. It consists of superconductor electrode 
(S), insulator (I) and FS bilayer as an upper electrode. 
We assume that the F film has a thickness, dp, and 
that it subdivides into domain structure with antiparal¬ 
lel direction of magnetization vector in the neighboring 
domains. The width of the domains is W and they sep¬ 
arated by atomically sharp domain walls oriented per- 
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pendicular to SF interfaces. Due to periodicity of the 
structures we, without any loss of generality, can per¬ 
form our analysis within its half of the period, that is 
from —W 12 to VF/2. This element is enlarged in Fig.l. 
It consists of two halves of domains and domain wall 
separating them. 

We will suppose that the condition of dirty limit is 
fulfilled for all metals and that effective electron-phonon 
coupling constant is zero in F material. We will assume 
further that either temperature T is close to the crit¬ 
ical temperature of superconducting electrodes Tq or 
the suppression parameters yss = RbS’^bn/Pf^f at 
SF interface is large enough to permit the use of the 
linearized Usadel equations in F film of the structure. 
We will characterize the FF interface (domain wall) by 
the suppression parameter 7 = 1 , and the suppression 
parameter ^bf = RbfAbf/Pf^f, which can take any 
value. Here Rbs,Rbf and Abn,Abf are the resis¬ 
tances and areas of the SF and FF interfaces, ^s, and 

= {Df/‘ 2 ttTcYA are the decay lengths of S, F ma¬ 
terials, while PS and pp are their resistivities. Dp is 
diffusion coefficient in the F metal. 

Under the above conditions the proximity problem 
in the SF part of SIFS junction (0 < a; < dp) reduces to 
solution of the set of linearized Usadel equations m-m, 

m 


f 52 52 1 ~ w 

f (9^ 52 1 ~ W 


where D = oj/pTc,G.± = |D| ±z/isgn(a;), h = H/nTc, 
H, is exchange energy of ferromagnetic material, w = 
7 rT( 2 n -|- 1) are Matsubara frequencies. The spatial co¬ 
ordinates in o, @ are normalized on decay length ^p. 
To write these equations we have chosen the x and y 
axis in the directions perpendicular and parallel to the 
SF plane and put the origin in the middle of SF inter¬ 
face to the point, which belongs to the domain wall (see 
Fig.l). 

Equations o, m must be supplemented by the 
boundary conditions |15) . They have the form 


d ^ A W W 

iBS^Fp = -Go — , a; = 0, - — <?/<—, 
ax oj 2 2 

d , W , X 

-Fp = 0, x = dp, --<y<- (3) 


At FF interface (y = 0, 0 < x < dp) and in the middle 
of the domains (y = ±Wl2, 0 < a; < dp) we also have 


1bf^Ff{x, -I-O) = Ff{x, -I-O) — Ff{x, —0), (4) 

dy 


■^Ff{x,+0) = ^Ff(x,-0), 

d W d W 

= ) = 0. 


( 5 ) 


Here W is the width of the domains, Gq = -I- A^, 

A is the modulus of the order parameter of supercon¬ 
ducting electrodes. The critical current density, Jc, of 
SIFS Josephson junction is determined by s-wave super¬ 
conducting correlations at IF interface, which is even 
function of the Matsubara frequencies 


sJcRn 

2nTc 


T 

WTc 


E 

CD>0 


GoA 




( 6 ) 


where $(y) = (FF,+a;(dF,y) -f Fp-(.,;(dF,y))/ 2 , while 
the full critical current, Iq, is the result of integration 
of Jciy) over width of the junction. 


cIcRn 

2nTc 


T 

WTc 



4>(y)dy. 


( 7 ) 


Here, Rn, is the normal junction resistance. 

Solution of Usadel equations in FS electrode. 
Solution of two-dimensional boundary value problem 
O-® in the F layer (0 < x < dp) is convenient to 
find in the form of the Fourier series expansion 

00 TT/’ 

FFix,y)= ^n(2/)cos-^, 0 < y < y, (8) 

n— — oo ^ 

Ff= ^ Bn{y) cos - Y < ?/ < 0, (9) 

n— — oo 

where 

Z / W\ 

An{y) =-^ + an cosh{q+ \^y - Y j), (10) 

Z / W\ 

Bn{y) = Y + bn cosh{q_ \y+ yJ), (H) 

and coefficients a„ and 




— 





iBsdpUJ 


are determined from boundary conditions Q. Here the 
coefficients <5, C± and S± are defined by expressions 


S = q-q+jBFS+S- -I- q-C+S- + q+S+C-, (14) 


G± = cosh(^^), Fi = sinh(^^). 


(15) 
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Taking into account the symmetry relation q-[—ui) = 
q+{^) for s-wave superconducting component in the F 
layer at x = it is easy to get 


<i>(2/>o) = - ^ (-ir 


n— — oo 


$(y< 0 ) = - ^ (-ir 


n— — oo 


■ I 

I 

■ I 

I ■ 

5+' 

E 

"(-2- 

q- 



5 _ 





(16) 

■ I 

I 

■ I 

I ■ 

<5_' 


"(-2- 

q- 

.E' 

~e. 

T 


(17) 

<3± = 9-J- cosh(g+-^-) - q+b+ cosh(g_- - -). 

Finally for the critical current from ©, (HU) and 
m we have 

cIcRn T 


zGqA 


2ttTc 2WTc -- w 


(18) 


S{uj) 


E(-i) 



W 


2S.S+(q^_ - qlf 
Sq^ q'j_ 


It is seen that the critical current can be represented 
as the sum of two terms. The first is the contributions 
from individual domains separated by fully opaque FF 
wall 


cIciRn _ T GqA^ 
2ttTc Tc ^ 

CJ .>U 


while the second 



eIc2RN ^ GgA" ^ {-1)^S-S+ 

2nTc WdpTc q+Q-S 

( 20 ) 

gives the contribution from the domain wall. Here Re(a) 
denotes the real part of a. 

Expression (HU) reproduces the well-known result 
previously obtained for single-domain SIFS structures 
m-m thereby demonstrating the independence of the 
critical current on the orientation of the domains mag¬ 
netization vectors, if they are collinear oriented and the 
FF interface is fully opaque for electrons. 

Limit of large jbf - For large values of suppression 
parameter 'jbf ^ max {l, (kFg±)“^} expression (1^ 
transforms to 


eIc2RN ^ GgA" ^ (-1)" 

2ttTc WdpTc ^ iBFlBS^'^ n^oo ^+®- 

( 21 ) 

The sum over n in Eq. (ED) can be calculated analyti¬ 
cally using the theory of residues 


cIc2Rn _ 2hT GgA^ 

2TrTc WTc ^ '^BFlBS^^ 

U) j>\J 


( 22 ) 


Si = Re 



/ 


1 


cosh 




-k 


dpy 


sinh [ dp y r2-|_ 


yj 


It is seen that Ic 2 is vanished as {''^bfW)~^ with 
increase of ^bfW product and scales on the same char¬ 
acteristic lengths ^ 1 , ^2 as the critical current for single¬ 
domain SIFS structures (ITOl) . 

Limit of small ^bf- In the opposite limit, ^bf "C 
max {l, (IF( 7 ±)“^} we have 


eIc2RN ^ 8h^T y. GgA^ 
27rTc WdpTc 


^ ^ {-irs.s+ 

" ^£^^qU^.{q-CpS.+q+S+C.y 

It is seen that in full agreement with the result ob¬ 
tained in m in the considered limit of large domain 
width, W Re(g±), 


eIc2RN GgA" ^ (-1)" 

2fTc WdpTc 7bsw2 q%ql_ (g_ -k q+)' 

(24) 

contribution to the critical current from domain wall 
region falls as W~^ and decays in the scale of fi. 

Limit of small domain width. In the opposite 
case, W <C Re(( 7 ±), presentation of the critical current 
as a sum of Ici and Ic 2 is not physically reasonable and 
for Ic from (1151) we get 


elcRjv GgA^ 

27rrq 2Tc jBsdFi^'^ 

to j>l} 


(25) 


Ss 


E (-1) 


(ql + g+) 7BW + 4 

(qlqllBW + q- + ql) 


where ysw = 7 b_fIF/2 . It is seen that for jbw ^ 1 ex¬ 
pression (ESI) transforms to (fTO|) and Ic = Ici, while in 
the limit ^bw 0 from (1251) it follows that the critical 
current 


el^ = il V-V (26) 

2TrTc Tc sinh 


is independent on exchange energy and falls with in¬ 
crease of dp in the same scale as it is for SINS devices. 
Previously it was found that such transformation of de¬ 
cay length takes place in a vicinity of domain wall |20] 
- m- In particular, it was shown that if a sharp do¬ 
main wall is parallel [52], [23] or perpendicular to SF 
interface [33] and the thickness of ferromagnetic lay¬ 
ers, df < ^p, then for antiparallel direction of mag¬ 
netization the exchange field effectively averages out. 
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and the decay length of superconducting correlations 
becomes close to that of a single nonmagnetic N metal 
= sjD p/2ttTc. The same effect may also take place 
in S-FNF-S variable thickness bridges [53], [35] . 

For arbitrary values of ysw the sum over n in (l25ll 
can be also calculated analytically. The denominator in 
(ESI) has the poles at 


TT y JEW 

Application of the residue theorem to the summation of 
the series in n in the expression (1551) leads to 


cIcRn 

2ttTc 


T "fBM ^ 


(27) 


^_ Q _ P _ 

^ yrr+psinh (di?-yrr+p) -yn + gslnh [dp^^l + q) ’ 

P ^ 1 - ^ ^ 1 +^ 28 ) 

IBW IBW 

It is seen that for "fBwh < 1 s-wave superconducting 
correlations decay exponentially into the F metal with¬ 
out any oscillations with two characteristic scales, = 
^F{n+p)~^/'^, and, ^12 = If 7BW tends 

to zero then one of the damping characteristic scale 
goes to that of SINF junctions (see (1551) 1. while 

the other fi 2 goes to zero. With jbw increase re¬ 
duces, whereas ^12 increases, so that at ypwh = 1 they 
become equal to each other = ^12 = ^f . 

Further increase of yBwh leads to appearance of the 
damped oscillations in Icidp) dependence with the ra¬ 
tio 

^ _ VIbw^'^ ~ ^ _ 

y/ ilBW^ + 1)^ + Ibw^"^ ~ 1 + ^IBW + 1 

(29) 

which monotonically increase from zero at yBwh = 1 
up to that of single domain SIFS junctions 


calculate the dependence of Icidp) is necessary to use 
the general expression dni). Figure 2 gives the Icidp) 
curves calculated for H = IOttFc, 7 bf = 0 and for a 
set of widths W/^p- It is seen that in full accordance 
with the analytical analysis given above for W smaller 
than 0.78 ^f, Ic falls monotonically with IF increase. At 
IF > 0.78 there is a transformation from a monotonic 
dependence of Icidp) to a damped oscillatory one. It is 
interesting to note that in the vicinity of the transition 
the critical current decays even faster than for large IF. 
To illustrate this result, we make a fit of the calculated 



FiglD Dependence of the critical current of SIFS 
Josephson junction as a function of thickness of F layer 
dp calculated numerically from (1181) for T = O.STc, 
H — IOttTc, 7sf = 0 and for a set of widths W/^f = 
0,3;0,5;0,7;0,8;1; 1,2. 


curves by the simple expression 

Icidp) = Aexp(-dF/6)cos(dF/C2 + p), 

which is ordinary used for estimation of the decay 
lengths and ^2 from an experimental data m. m- 
At the first step we define ^2 


in the limit Jew 00 . 

From ESI) ESI) we can conclude that the existence 
of domain structure in the F layer of SIFS devices can 
significantly modify the relation between ^1 and ^2 ex¬ 
tracted from experimental studies of Icidp) dependence 
in SIFS tunnel junctions. 

This conclusion is valid not only in the limit of small 
domain width. 

Arbitrary values of the domain width. For ar¬ 
bitrary values of the width of the magnetic domains to 


^2 = idF2 — dFi)/Tr 


from the positions of the first, dpi, and the second, dF 2 , 
O-TT transitions in Icidp) dependence and put 


P = 7r/2 - dFi/6 


in order to get Icidpi) = 0. The decay length is de¬ 
termined from the ratio of magnitudes of critical current 
taken in two points having equal phase of oscillation: 


Cl = 7^6 In 


Icidpi + C27r/2) 
IcidF2 + C27r/2) 

































Josephson effect in SIFS tunnel junctions with domain walls in weak link region 


5 


and normalization constant A 

^ _ Icidpi + 67r/2) _ 

exp(-dF/Ci) cos(dF/6 + ’f) 

has been determined by direct calculation of magnitude 
in the certain point between O-tt transitions. If the po¬ 
sition of the second O-tt transition exceeds 10 we 
suppose that ^2 is infinite and Ic{dFi) dependence can 
be fitted by function 

Jc{dF) = Aeyip{-dFlii). 

The results of the fitting procedure are presented in 
Fig.3-Fig.5, which give the decay lengths and ^2 
as well as their ratio Ci /?2 calculated at T = O.STc, 
F[ = IOfTc for a set of suppression parameter jff = 
0;0.3;1. Thin vertical lines in Fig.3, Fig.4 give values 
on the x-axis, at which there is a transition from a 
monotonous exponential decay of Ic[dF) to the damped 
oscillation lows. Thin horizontal lines in Fig.3 - Fig.5 
provide the asymptotic values of ^ 1 , ^2 and ^ 1/^2 in 
the limit VF which are coincide with the mag¬ 

nitudes calculated for single domain SIFS junction for 
given temperature T = 0.5Tc and exchange energy 
H = IOttTc. 



FigEl Dependence of decay length as a function of 
domain width W calculated at T = O.STc, H = IOttTc 
and 7 sf = 0; 0.3; 1. 


It is seen that the transition point at which mono¬ 
tonic decay of Ic(dFi) dependence transforms to a 
damped oscillation behavior the smaller the larger is 
suppression parameter 'ysF- Interestingly, in the vicin¬ 
ity of this transition decay length is even smaller com¬ 
pare to its magnitude in the limit of large W. 

It is also necessary to note that despite of the fact 
that the transition takes place at IF < ^f, the differ¬ 
ence between and ^ 2 , as it follows from Fig.5, exists 



FigEl Dependence of decay length ^2 as a function of 
domain width IF calculated at T = O.STc, H = IOttTc 
and 7 sf = 0; 0.3; 1. 



tion of domain width IF calculated at T = O.STc, 

H — IOttTc and 7 bf = 0; 0.3; 1. 

even for large domain width: the ratio C 1 /C 2 is only 
around 0.8 at IF = 4,^^ and very slowly tends to the 
following from (15(11) the single domain value 0.95 with 
IF increase. This fact permits us to conclude that the 
difference between and ^2 experimentally observed in 
SFS Josephson structures based on dilute magnetic al¬ 
loys can be also the consequence of existence of magnetic 
domains in the F layer. 
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Fig.l. Geometry of the considered SIFS Josephson 
junction and its enlarged part, which includes two halves 
of domains and domain wall separating them. The in¬ 
sulating barrier I has a small transparency (shown by a 
blue line). 

Fig.2. Dependence of the critical current of SIFS 
Josephson junction as a function of thickness of F layer 
dp calculated numerically from (TT^ for T = 0.5Tc, 
H = IOttTc, 7 bf = 0 and for a set of widths W/^p = 
0,3;0,5;0,7;0,8;I;1,2. 

Fig.3. Dependence of decay length as a function of 
domain width W calculated at T = 0.5Tc, H = IOpTc 
and ^bf =0; 0.3; I. 

Fig.4. Dependence of decay length ^2 as a function of 
domain width W calculated at T = 0.5Tc, H = IOttTc 
and 7 bf =0; 0.3; I. 

Fig.5. The ratio of decay lengths and ^2 as a 
function of domain width W calculated at T = 0.5Tc, 
H = IOttTc and ')bf = 0; 0.3; I. 



